Entanglement of Fermi gases in a harmonic trap 
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For both cases with and without interactions, bipartite entanglement of two-fermions from a Fermi 
gas in a trap is investigated. We show how the entanglement depends on the locations of the two 
fermions and the total number of particles. Fermions at the verge of trap have longer entanglement 
distance (beyond it, the entanglement disappears) than those in the center. We derive a lower 
limitation to the average overlapping for two entangled fermions in the BCS ground state, it is 
shown to be yjQ/2M, a function of Cooper pair number Q and total number of occupied energy 
levels M. 
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As one of the most characteristic features of quan- 
tum systems, quantum entanglement lies at the heart of 
the difference between the quantum and classical multi- 
particle world. It is the phenomenon that enables quan- 
tum information processing and quantum computing 0, 
2] . Quantum entanglement is usually considered as exist- 
ing between different degrees of freedom of two or more 
particles with mutual interactions, it is only recently that 
researchers have started to investigate entanglement in 
systems containing a large number of particles, in par- 
ticular in a noninteracting Fermi gas Q-OI-Q- Entangle- 
ment seems to play a crucial role in condensed matter 
systems, and has shown of relevance to thermodynami- 
cal quantities such as the degeneracy pressure|6| and the 
number density of the gas|5j, multipartite entanglement 
is also promising to make breakthrough in solving un- 
solved problems such as high-T c superconductivity Q- 

Entanglement for noninteracting Fermi gases in a free 
space has already been studied0,0,0|- It was found that 
all entanglement vanishes if the relative distance \r — 
r'\ between electrons is greater than the entanglement 
distance 1/kp, where fop denotes the Fermi momentum. 
In this situation, quantum entanglement is purely due to 
Fermi statistics and not due to any physical interaction. 
A natural question then arises, is this a general property 
of noninteracting Fermi gases? or it is special for the 
Fermi gas in a free space, and how the interactions among 
the fermions influence the entanglement? 

In this paper, we will try to answer this question by 
studying the entanglement in noninteracting Fermi gases 
trapped in a harmonic trap, and examining the effect of 
inter- fermion interactions on the entanglement. As you 
will see, bipartite entanglement measured by Wootters 
concurrence depends not only on the relative distance 
between two fermions but also on the total number of 
particles, the larger the number of fermion N, the shorter 
the entanglement distance. The entanglement distance is 
no longer a constant as in the free space, it is related 
to locations of the two fermions. Our numerical simula- 
tions show that the entanglement distance is longer at the 
verge more than that at the center of the trap. Further, 
we show the effect of interactions on the bipartite entan- 
glement for the Fermi gas. The model adopted is the 



reduced BCS Hamiltonian, it is shown that the bipartite 
entanglement of BCS ground state sharply depends on 
the overlapping of the time-reversed states, the minimum 
average overlapping is y / Q/2M depending on the ratio 
of Cooper pair number Q to the total number of occupied 
energy levels M, this indicates bipartite entanglement in 
the conventional BCS ground state. These results clearly 
establish the fact that entanglement should be taken into 
account when studying macroscopic observable, even if 
the system is in its ground state and the constituents of 
the system are noninteracting as the fermions in a trap. 

Suppose we have a collection of noninteracting Fermi 
gas in harmonic traps. The Fermi gas may be electrons 
in a metal or ultracold atoms. Unless stated otherwise, 
we treat the Fermi gas in this paper as ultracold atoms, 
but the representation is applicable for all Fermi gases. 
At zero temperature, the atom gas is in its lowest en- 
ergy configuration. All energy states are occupied up 
to the level M = N/2, where N is the total number of 
atom and assumed to be even (the case with odd number 
of particles will be discussed at the end of this paper). 
For simplicity, we consider throughout this paper a one- 
dimensional harmonic trap. This condition is met if the 
trapping frequencies in the other directions are consider- 
able large. The ground state of this system is 

M 

l*o> = Wbl^Vac), (1) 

n=l 

where \Vac) denotes the vacuum, and b} l(J (a =t,|) cre- 
ates an atom in state <fi n (%) with spin a. We are inter- 
ested in entanglement between two atom spins at differ- 
ent locations. The solution to this problem would answer 
the following question, suppose that one atom has spin 
up at x, can we infer the direction of the spin of the other 
atom at x'l In order to answer this question, the density 
matrix describing the spin state of two atoms at locations 
x and x' is needed. It can be defined by§ 

P..>,tt> = (* |4(^(s)Mz>.(z)|*o>, (2) 

where ip}(x) creates an atom of spin t at the location x. 
This density matrix also may be calculated by 

p ss ,,u> = Tr(|*o)<*o| • \st(x)){s't'(x')\), (3) 
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with \st(r)) standing for two-atom state with spins s and 
t at the location r. Writing ij) s (t) in terms of the anni- 
hilation operator b na and the eigenfunctions of harmonic 
oscillator (j) n (x), we obtain the density matrix in the fol- 
lowing form, 

Pss ,, tv = N(x)N(x')5 ts 5 t > s > - S ts ,5 t , s F 2 (x, x'), (4) 



where F(x,x') = E"ftMfcW = F, N(x) = 
^2 a \<fi a (x)\ 2 — N x . Term F represents a sum over over- 
lapping of the two atoms at locations x and x' , respec- 
tively. N(x) is the number density of atom at the location 
x. In basis {| ft), I Tl), I IT), I II)}, the density matrix 
takes the form, 



pu(x,x') = 
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FIG. 1: Wootters concurrence of two atoms located at x and 
x' , respectively. The figure was plotted for 20 trapped atoms. 
a — -J mu/h, m is the mass of atom, and lo is the trapping 
frequency. 



FIG. 2: The bipartite entanglement as a function of the sec- 
ond atom's location x. (a)-(e) are for different total num- 
ber of atoms. The number corresponding to (a),...,(e) is 
2,4,10,14,18, respectively. 



where the subscript 1,2 indicates that there are two 
atoms. The entanglement of formation [j| measured by 
Wootters concurrence can be given by 

Ci2(x,x') = mxN l_ 2F 2\ ™*{ 2F2 - N * N *'^}- ( fi ) 

Obviously, the Wootters concurrence C12 is maximal 
when x = x' and it equals to 1. The corresponding en- 
tanglement state is the spin singlet l/-\/2(| ||) — | IT))- 
With the relative distance between the two atoms in- 
creasing, F 2 behaves as a damping function of \x — x'\. 
As a consequence, the entanglement decay with \x — x'\ 
increasing, this was shown in figure 1, where the bipartite 
entanglement in a system with 20 atoms was plotted as 
a function of locations x and x' of the two atoms. The 
entanglement arrives at the maximum value 1 at points 
I a; — x'\ — and decays with atom separation increasing 
from the pauli exclusion principle. This can be under- 
stood as a result of more and more triplet states mixed 
in with the singlet. In contrast with the entanglement 
between two fermions in a free space, the entanglement 



between the two in a harmonic trap appears to be lo- 
cation dependent. It is clear from figure 1 that the en- 
tanglement distance is longer for atoms at the verge of 
the trap than that at the center. The reason for this is 
the following. The bipartite entanglement roughly de- 
pends on how and how many triplet state are mixed in 
with the singlet. At points \x — x'\ — 0, there are no 
triplet states involved in from the pauli exclusion, so the 
bipartite entanglement is maximal at these points. With 
the atom separation increasing, more triplet states are 
involved in, the two atom state is then a weighted sum 
over the singlet and triplet states, the weights depend on 
the locations of atoms. Because the center of the trap is 
the favorite location for particles to occupy, the weights 
at the verge benefit the entanglement. The bipartite en- 
tanglement also depends on the total number of atom, as 
shown in figure 2, where one atom was located at 0.5a, 
a = \J muj/h, m is the mass of atom, and u> stands for 
the trapping frequency. Clearly, the larger the number 
of atom, the faster the damping of the bipartite entan- 
glement. 
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Up to now, we have not considered the interaction 
between the atoms. In the following, we study how 
the interaction influences the entanglement of the two 
atoms. For repulsive interactions, perturbation theory 
tells us that weak interactions modify the eigenvalues and 
the corresponding cigcnfunctions of the free Hamiltonian, 
this results in a level shift to every eigenstate of the free 
Hamiltonian, and consequently there are more eigenfunc- 
tions participating in the summations in F, N x , and N x >- 
So, the damping in the dependence of entanglement on 
the relative distance would faster than that without in- 
teractions. For attractive interaction, we will investigate 
the entanglement in trpped atoms by using the reduced 
BCS model. The reduced BCS Hamiltonian has received 
much attention as a result of effort to understand pair- 
ing correlations in nanoscale metallic system |9j, Il0| . The 
model Hamiltonian is[TT| . 

M M 

Hbcs = - d\^b] + b)_b k+ b k -, (7) 

where j, k = 1,...,M labels a set of doubly degenerate 
single particle energy levels with energies ej, A is the 
dimensionless coupling, and d is the mean level spac- 
ing. bj±(bj±) represent the annihilation (creation) op- 
erators for atoms at level j with the labels ± referring 
to pairs of time-reversed states, and nj was defined as 
rij = b J j + bj + + bj_bj-, the atom number operator for 
level j. The ground state entanglement (called ALC - 
average local concurrence) share among different energy 
levels in this model was investigated in (lj), it shown 
a simple relation between the ALC and the order pa- 
rameter. The conventional BCS theory employs a grand 
canonical ensemble, by the Bogoliubov transformation 



7ji = Ujb j+ -Vjb!j_, 7jo = Ujbj- +vjbj + , with u 2 + v 2 

I. 4u}tf = A 2 /(e 2 + A 2 ), and A = \dJ2*L 1 {bj-b j+ ), it 
gives 

Hbcs = Yl^^o ~ v ^ ~ 2An i^K7ji7ji + 7jo7jo) 
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+ constant. (8) 

Instead considering entanglement between two atom 
spins, we here consider entanglement between two atoms 
in the two time-reversed states + and — . The density 
matrix representing this entangled state can be defined 
as 

pf s ? : ? t , = Tr(\BCS){BCS\ ■ \st(x))(s't' (x')\), (9) 

where \BCS) is the well-known ground state in the BCS 
model, \BCS) = Ylf^M + v 3 b] + b]^)\Q), s,t,s',t' = 
+, — . By the standard precedure, Eq. © yields, 

PS = (Y,h\ 2 ) 2 S ts 8 t , s/ -S ts ,S t/s Re(f(x, X iy(x,xi)), 

(10) 

where 

f(x,x') = ^^(x)0*(^') = / 5 
3 

v 2 (x,x') = J2 v ]<l>*i( x ) ( t>j(x') = v 2 . (11) 
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Here, Re(...) denotes the real part of (...). In basis {| + 

+), | H — ), | — h), | )}, the density matrix defined in 

Eq.@ follows, 
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(12) 



212 
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v 2 equals the total number Q of Cooper pairs in 
state \BCS), while / characterize the overlapping of the 
two particles. The state Eq. l|12|l is entangled iff the 
Peres-Horodecki condition is met, i.e., 2Re(fv 2 ) — 
Q 2 > 0. Assume (f)*(x)4>j(x') — y, y is a constant, 
the Peres-Horodecki condition leads to |y| > y/Q/2M 
(< l/y/2). It is the restriction to the overage overlap- 
ping for bipartite entanglement. With the assumption 
<fij(x)<f)j(x') = y, it is easy to write down the concurrence 

Ci 2 CS {x, x') = max{ [^q/m-i^P) » Q}- The maximal entan- 
glement C^ 2 cs (x,x') = 1 is obtained at \y\ 2 = Q/M. 



For conventional BCS state, </>j(x) may have the form of 
e lp i x , choose x = x' , \y\ = 1 > l/v2- So two electrons 
in the conventional BCS state are entangled as long as 
the separation L is less than the enta ngleme nt distance 
satisfying \4>*(x)4>j(x + L)\ = \y\ = y/Q/2M. However, 
it is unsure for a system with (f>i{x) taking the eigen- 
states of harmonic oscillator. It is easy to show that 
\<p*(x)(j)j(x')\ < \j\[2 for any j, x and x', which indi- 
cates that small ratio of the Cooper pair number to the 
total number of occupied levels might guarantee the en- 
tanglement existing in such a system. 

We now consider the case when the total number of 
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particles is odd. Intuitively, for a large number (say, 27V + 
1) of atom, the bipartite entanglement would behave like 
that with 27V particles. This is the case indeed as you will 
see. Let us first analyze the case without inter-particle 
interactions. Assume the (27V + l)th atom is of spin 
up, the state of this 27V + 1 atoms may be written as 
b\ I+1 ^o), where |W ) was defined by Eq.JTJ for the 27V 
particles. Following the calculation performed for Eq.(J2J), 
we find that the density matrix pi2(x, x') can be divided 
into two matrices, the first represents contributions from 
the 27V atoms, which takes the same form as in Eq.©, 
and the second is a correction due to the (27V+l)th atom. 
The elements of the second matrix a are, 

a 22 = 2N x \(f> M+ i(x')\ 2 + 2N x ,\(f> M+1 (x)\ 2 , 

<T23 = -4F(x,x')<t>* M+1 (x)4> M +l(x'), (13) 

oil = tf22 +CT23, 022 = 0-33, 032 = 0"f 3 , and the others are 
zero. For a large system (TV >> 1), \(j)M+i{x)\ 2 « N x , 
and I^m+i^^az+i^OI << \F{x,x')\. Therefore, the 
correction a to the density matrix p\2(x,x') can be ne- 
glected. However, this is not the case if the system only 
consists of few atoms, the correction due to the (27V+l)th 
atom should be taken into account to compute the bipar- 
tite entanglement. Note that exchanges of spin up with 
spin down do not change l^o); but fc] v/+1 ||^E'o)- So ma- 
trix a is of relevance to the spin of the (27V + l)th atom. 
But this does not affect the bipartite entanglement under 
consideration, i.e., the bipartite entanglement is indepen- 
dent of the spin of the (27V + l)th atom. In the case of 
attractive interaction, the situation is similar if the BCS 
ground state is simply b\ [+l ^\BCS). The situation be- 
comes complicated when b\ I+l ABCS) is not the ground 



state of the system 14] . In order to calculate the entan- 
glement in BCS ground state, we have to find the ground 
state first, it is beyond the scope of this paper. 

Before concluding, it is worth mentioning that bi- 
partite (multipartite) entanglement as a properties be- 
tween quantum systems depends on definition of the two 
degrees (many degrees) which share the entanglement. 
Some properties of entanglement in the BCS model were 
studied in Ref. 0, 0] , where the entanglement was de- 
fined among particles in different energy levels ^3 ■ The 
entanglement presented in this paper is for two fermion 
spins at different locations x and x' , it characterizes two- 
fermion's correlation at the two locations. 

In conclusion, we have shown that the bipartite en- 
tanglement in noninteracting fermions trapped in har- 
monic traps depend on particle number, relative distance 
and the locations of the two fermions. The entangle- 
ment distance which characterize the maximum separa- 
tion of two entangled particles is longer at the verge of 
the trap than that at the center; the larger the number 
of trapped particles, the shorter the entanglement dis- 
tance. For interacting Fermi system, we have adopted 
the reduced BCS model to study the entanglement in 
the BCS ground state. Reduced density matrix and 
Wootters concurrence have been presented, the restric- 
tion on the average overlapping y has been derived to be 
y/Q/2M < \y\ < y/Q/M, the lower limitation y/Q/2M 
co rrespo nds to concurrence zero, while the upper one 
y/Q/M corresponds to maximal entangled states. 
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